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ABSTRACT 

We use models of thermal evolution and XUV-driven mass loss to explore the composition and history of 
low-mass low-density transiting planets. We investigate the Kepler- 11 system in detail and provide estimates 
of both the current and past planetary compositions. We find that a H/He envelope on Kepler- lib is highly 
vulnerable to mass loss. By comparing to formation models, we show that in situ formation of the system 
is extremely difficult. Instead we propose that it is a water-rich system of sub-Neptunes that migrated from 
beyond the snow line. For the broader population of observed planets, we show that there is a threshold in bulk 
planet density and incident flux above which no low-mass transiting planets have been observed. We suggest 
that this threshold is due to the instability of H/He envelopes to XUV-driven mass loss. Importantly, we find 
that this mass loss threshold is well reproduced by our thermal evolution/contraction models that incorporate 
a standard mass loss prescription. Treating the planets' contraction history is essential because the planets 
have significantly larger radii during the early era of high XUV fluxes. Over time low mass planets with 
H/He envelopes can be transformed into water-dominated worlds with steam envelopes or rocky super-Earths. 
Finally, we use this threshold to provide likely minimum masses and radial velocity amplitudes for the general 
population of Kepler candidates. Likewise, we use this threshold to provide constraints on the maximum radii 
of low-mass planets found by radial velocity surveys. 

Subject headings: planetary systems; planets and satellites: composition, formation, interiors, physical evolu- 
tion; stains: individual (Kepler- 11) 



1. INTRODUCTION 

In recent years, the frontier of the search for extrasolar 
planets has pushed towards ever smaller and more Earth-like 
worlds. We now know of dozens of Neptune mass planets 
and hav e even found the first definitively ro cky extrasolar 
planets dBatalha et all 1201 U ITeger et al.ll2009h . In between, 
transit searches have begun finding a population of low-mass 
low-density "super-Earths". Beginning with the discovery of 
GJ1214b (Charbonneau et al. 2009), these planets represent a 
new class of exoplanets that do not have any analog in our 
Solar System. Basic questions about their composition, struc- 
ture, and formation are still unknown. Are these, in fact, 
scaled up versions of the Earth that simply have thick hydro- 
gen/helium envelopes atop of rock/iron cores? Or are they 
instead scaled down versions of Neptune that are rich in wa- 
ter and other volatile ices? 

The distinction between water-poor super-Earths or water- 
rich sub-Neptunes has fundamental implications for how 
these planets formed. So far these low-mass low-density 
(hereafter LMLD) planets have only been found well inside 
the snow-line. If these planets only contain rock, iron, and 
hydrogen/heli um, then it is possible the y formed close to their 
current orbits ([Hansen & Murravl201 ll) . However, if a signif- 
icant fraction of their mass is in water, then they must have 
formed be yond the snow-line and migrated in to their current 
locati ons (lAlibert et al.ll201 U llda & Linll20Tot [Rogers et al.1 
l20Tlh . 

The Kepler- 11 system jLissauer et al.l |201 la) is an ex- 
tremely powerful tool for exploring the features of LMLD 
planets. With six transiting planets orbiting a close solar ana- 
log, it is the richest extrasolar system currently known. More- 
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over, five of the planets have masses from Transit Timing 
Variations (TTVs), and all five of these fall into the low-mass 
low-density regime in between Earth and Neptune. These 
five planets are all interior to Mercury's orbit, with periods 
from 10 to 47 days. This provides a unique laboratory to test 
the possible composition, formation, and evolution of LMLD 
planets and how these vary as a function of both period and 
planet mass. 

Transiting planets with measured masses, like those in 
Kepler- 11, are particularly valuable because we can determine 
their mean density. All the planets in Kepler- 1 1 have densities 
too low for pure rock, and therefore must have some sort of 
thick envelope of volatiles. Likewise, all the planets except 
Kepler- 1 lb are less dense than pure water and so must have at 
least some hydrogen/helium. 

Unfortunately, mass and radius alone cannot uniquely de- 
termine a planet's composition. In general, there is a large 
degeneracy between th e relative amounts of rock , iron, water, 
and hydrogen/helium ([Rogers & Seagerll2010ab . This prob- 
lem is particularly acute for planets with radii « 
since in this range any of these four constituents can be im- 
portant. Indeed these sorts of degenerac ies have long been a 
focus o f studies of Uranus and Neptune ([Hubbard et al.[[T99Tt 
[Fortnev et al. 2011). 

One possible solution to the composition problem is to ob- 
tain multi- wa ve length transmission spectra, as has been done 
for GJ1214b ([Bean et al.l[20TTt [Desert et al.[[201 U [Croll et al.1 
l20Tlh . Since hydrogen-rich atmospheres have much larger 
scale heights at a given temperature, near infrared water and 
methane absorption features will be much mo re prominent for 
planet s with hydrogen/helium env elopes ( Miller- Ricci et al.1 
2009: Miller-Ricci & Fortnevl2010l) . Unfortunately, these ob- 
servations are extremely time intensive and even then the pos- 
sible presence of clouds can make their interpretation diffi- 
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cult. Even worse, nearly all the systems found by Kepler are 
too faint for these observations with current telescopes. 

An alternative is to develop models of the formation and 
evolution of low-mass planets to try and predict what com- 
positions can form and how those compositions change as a 
planet evolves. In particular, hydrodynamic mass loss from 
extreme ultra-violet (XUV) heating can remove large amounts 
of hydrogen/helium from highly irradiated LMLD planets. 
Models of XUV driven mass loss were first developed to study 
water loss from early Venus ('H untenll982l lKasting & Pollack' 
1983 ), and h ydrogen loss from the early Earth (Sekiva et al. 
T98a Imtson et all IT981I) . These kinds of models have 



since been developed to st udy m ass loss from hot Jupiters 
(e.g., p^ammer et al. 2003^ 
l2009t lEhrenreich & Dese^ 
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12004; Murrav-Clav et al.] 
Owen & Jackson 2012), 
where there is strong evidence that atmospheric escape 
is an im portant physical process (Vidal-Madjaret al. 2004; 
Davis & Wheatlev 2009; Leca velier Pes Etangs et al. ,2010; 
Lecavelier des Etangs e t al. 201 2|k 

In Sections 12.31 13.31 and 15.11 we show that energy-limited 
hydrodynamic mass loss models, coupled with models of 
thermal evolution and contraction, can distinguish between 
water-poor super-Earth and water-rich sub-Neptune scenarios 
in Kepler- 1 1 . Moreover, these models make powerful pre- 
dictions for the density distribution of the entire population 
of LMLD transiting planets. In particular, observations show 
that there is threshold in the bulk density - incident flux distri- 
bution above which there are no LMLD planets. In Section|4] 
we examine this threshold and show how it can by reproduced 
using our thermal evolution models coupled with standard hy- 
drodynamic mass loss prescriptions. Finally, in Section 14.21 
we explore how this threshold can be used to obtain impor- 
tant constraints on planets without measured densities: We 
constrain the maximum radii of non-transiting radial velocity 
planets, and the minimum masses of Kepler candidates. 

2. OUR MODEL 

2.1. Planet Structure 



We have built on previou s work in iFortney et al.l (l2007h 
and Nettelman n et al.l (1201 Ih to develop models of the ther- 
mal evolution of LMLD planets. To simplify what is undoubt- 
edly a complex interior structure for real planets, we construct 
model planets with well-defined layers. Low-mass planets are 
likely to have a significant fraction of their mass in iron and 
silicate rocks. For simplicity, we assume that these materials 
are contained in a isothermal rocky core with Earth-like pro- 
portions of 2/3 si hcate rock and 1 /3 iron. For the rock, we 
use the ANEOS dThompsonl 1 1 990h olivine equation of state 
(EOS); while for the iro n, we use the SESAME 2140 Fe EOS 
jLvon & JohnsonlfT992l) . 

On top of this rock/iron core we then attach an interior adi- 
abat. The composition of this adiabat depends on the planet 
model being considered. For this work, we consider three 
classes of LMLD planets: rocky super-Earths with H/He en- 
velopes, water-worlds that have pure water envelopes, and 
sub-Neptunes with a water layer in between the core and the 
upper H/He layer For the water-rich sub-Neptune models we 
assume that this intermediate water-layer has the same mass 
as the rock/iron core. We choose this value because it is 
comparable to the water to rock ratio need to fit Kepler- lib 
as a water-world. This allows us to explore the proposition 
that all five Kepler- 11 planets started out with similar com- 
positions, but that mass loss has subsequently distinguished 



them. For hydrogen/helium we use the ISaumon et al.l (11995b 
EOS. Meanwhile for water we use the ab-i nitio H20-RE0S 
EOS developed by 'Nettelmann et alj (l2008h and [French et all 
(2009), which was recently confirmed up to 7 Mbar in labora- 
tory experiments (Knudson et al. 20T2|). 

In the Kepler- 1 1 system, our models predict that water will 
be in the vapor, molecular fluid, and the ionic fluid phases. 
The interiors are too hot for high pressure ice phases. Finally, 
we model the radiative upper atmosphere by assuming that 
the planet becomes isothermal at pressures where the adiabat 
is cooler than the planet's equilibrium temperature, assuming 
20% Bond albedo and uniform re-radiation. We then calculate 
the radius at 10 mbar which we take to be the transiting radius. 

We connect the different layers of our models models by 
requiring that pressure and temperature are continuous across 
boundaries. We then solve for the interior structure assuming 
hydrostatic equilibrium. A given model is defined by its mass, 
composition (i.e., the relative proportions in H/He, water, and 
the rock/iron core) and the entropy of its interior H/He adia- 
bat. By tracking changes in composition and entropy we can 
then connect these models in time and study the thermal and 
structural evolution of a given planet. 

2.2. Thermal Evolution 

In order to obtain precise constraints on composition, it is 
important to fully model how a planet cools and contracts due 
to thermal evolution. Models that only compute an instanta- 
neous structure (Rogers & Seager 2010b) by necessity must 
vary the intrinsic luminosity of the planet over several or- 
ders of magnitude, which can introduce large uncertainties in 
the current composition. Obtaining precise constraints from 
thermal evolution is essential when considering mass loss, 
since mass loss histories are highly sensitive to uncertainties 
in the current composition. Moreover, since mass loss de- 
pends strongly on planetary radius (to the third power), the 
mass loss and thermal histories are inextricably linked. 

Modeling this contraction requires a detailed understanding 
of a planet's energy budget. By tracking the net luminosity of 
a planet, we know how the specific entropy S, (i.e., the en- 
tropy per unit mass) of the interior adiabat changes with time. 
For a given mass and composition, this adiabat then defines 
the planet's structure and so we can track the planet's total ra- 
dius as the model cools and contracts with time. Equation dU 
shows the energy budget for our models and how this relates 
to the change in entropy dS/dt. 



dm 



TdS 
dt 



— ^int "t" ^adio ^v^ct 



dTr-r 



dt 



(1) 



The left hand side shows the rate of change of the ther- 
mal energy of the interior adiabat. Positive terms on the right 
hand side represent energy sources that heat and inflate a 
planet, while negative terms represent energy losses that al- 
low a planet to cool and contract. The term Lint = i-eff - i-eq 
describes the intrinsic luminosity due to radiation from the 
planet, where Leq is the planet's luminosity due only to ab- 
sorbed stellar radiation. 

The Lradio term describes heating due to radioactive de- 
cay. The important isotopes are ^''^U, '^''K, ^^^U, and ^^^Th. 
These have half lives of 0.704, 1.27, 4.47, and 14.1 Gyr, 
respectively. We assume meteoritic abundances given by 
lAnders & Grevessel (119891) . We do not consider the early de- 
cay of ^^Al, since we only consider models that are at least 10 
Myr after planet formation. The Z^^adio term has only a minor 
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effect on our models since it is typically an order of magnitude 
smaller than the other terms in equation ([T). 

Lastly, there is the dTcoTc/dt term, which represents the de- 
lay in cooling due to the thermal inertia of the rocky core. As 
the interior adiabat cools, the core isotherm must also cool, 
as Tcore equals the temperature at the bottom of the adiabat. 
When the core makes up a large fraction of the planet's mass, 
this can significantly slow down the planet's rate of contrac- 
tion. We assume a core heat capacity of Cy = 0.5-1 .0 JK~' g~' 
dAlfe et al. 2002; Guillot et al.lll99l iValencia et all lMO*) as 
in ,Nettelmann et al. (20111). This range covers values appro- 
priate for both the cores of the Earth and Jupiter For our 
three layer sub-Neptune models we still us the mass in rock 
and iron for Mcok, since the water layer is generally too hot 
for ice phases and so it is assumed to be fully convective. 

For a given interior structure, we determine the intrinsic 
flux from the interior, at given S of the adiabat, via interpo- 
lation in a grid of model atmospheres. The values of Tim (a 
parametrization of the interior flux), Tgq, and Teff are tabulated 
on a grid of surface gravity, interior specific entropy, and inci- 
dent flux for 50 X solar metallicity H/He atmospheres (similar 
to Neptune). This corresponds to a metal mass fraction of 
Z « 0.35 and a mean molecular weight of /i 3.5 gmol"'. 
The g rid is the same as that described in iNetteknann et al.l 
( 1201 11) for LMLD planet GJ 1214b, where a more detailed 
description can be found. Here we do expand on that grid to 
now includ e a range of incident fluxes, as was done for giant 
planets in Fortne^et al. (2007). 

In choosing the initial entropy for our evolution model, we 
assume a "hot start" for model; i.e., we start the models out 
with a large initial entropy. When then allow the models to 
cool and contract until either 10 Myr or 100 Myr which is 
when we begin the coupled thermal and mass loss evolution. 
This is a common but important assumption. However, in 
general our thermal evolution models are insensitive to the ini - 
tial entropy choice by ^ 100 Myr as in iMarley et al.l (l2007l) . 
As a result, we present results at both 10 and 100 Myr More- 
over, to gain confidence in our 10 Myr models we examined 
the effect of starting those models with a lower initial entropy. 
Specifically, we ran models in which we started the 10 Myr 
with the entropies found at 100 Myr This allowed us to sep- 
arate the effect of the stellar XUV evolution, from any "hot 
start" vs. "cold start" uncertainties. Future progress in mod- 
eUng the formatio n of water-rich sub-Neptune planets (e.g., 
[Rogers etal.ll2011b may allow for an assessment of the most 
realistic initial specific entropies. 



2.3. XUV-Driven Mass Loss 

Close-in planets like those in Kepler- 1 1 are highly irradi- 
ated by extreme ultraviolet (EUV) and x-ray photons. These 
photons photoionize atomic hydrogen high in a pla net's at- 
mo sphere, which in turn produces significant heating (iHuntenI 
Il982[) . If this heating is large enough, it can generate 
a hydrodynamic wind that is capable of removing signif- 
icant mass, pot entially including heavier elements as well 
dKasting & Poll ack 1983). We couple this XUV driven mass 
loss to our thermal evolution models following the approach 
of Jackson e t al, ( 2010) and Valencia et al. (2010), wh ich ex- 
plored possible mass loss histories for CoRoT-7b (ILeger et alJ 
120091; lOueloz et al.ll2009l) . Similar approaches have also been 
used to study the coupled evolution of hot Jup iters (e.g., 
iBaraffe et a l.' "SOO?, '2005'; 'HubbardetaD l2007biah and hot 



A common approach to estimate the mass loss rate is to as- 
sume that some fixed fraction of the XUV energy incident on a 
planet is converted into heat that does work on the atmosphere 
to remov e mass. This is know n as the energy-limited approx- 
imation ("Watson et al.lll98 T) and allows a relatively simple 
analytic description of mass loss rates. 



Me. 
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/^tide-(l-2|+2^) 
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(3) 



(4) 



Neptunes fBaraffe et al. 2006.) . 



Equation dU describes our estimate of the mass l oss rate 
based on the formulation from (.Erkaev et al.l I2007h . Fxvw 
is the total flux between 1 - 1200 A, which is given by 
Ribas(2005) for Sun-like stars. For stars older than 100 Myr, 
Ribas found that at 1 AU fxuv = 29. 7t"''^-' ergs"' cm"^, where 
r is the age of the star in Gyr Using this power law, we scale 
the XUV flux to the appropriate age and semi-major axis for 
each planet in our models. Although Ribas only targeted Sun- 
like stars, Sanz-Forc ada et all (1201(1 ) found similar results for 
a wide range of stellar types from M3 to F7. Hereafter, we 
will simply refer to the entire 1 - 1200 A; spectrum as XUV. 

^xuv is flie planetary radius at which the atmo- 
sphere becomes optically thick to XUV photons, which 
iMurrav-CIav et all (120091) find occurs at pressures around 
a nanobar, in the hot Jupiter context. For our work, we 
assume that the atmosphere is isothermal between the op- 
tical and XUV photospheres. This neglects heating from 
photo-disassociation, which shou ld occur around a /xbar 
dMiller-Ricci Kempton et all 120121) . However, this effect 
should be relatively small and if anything will lead to slight 
underestimate of the mass loss rate. We vary pressure of the 
XUV photosphere from 0. 1 nbar to 10 nbar to include the un- 
certainty in the structure of the XUV photosphere. For H/He 
atmospheres on LMLD planets, the nbar radius is typically 
10-20% larger than the optical photosphere. /Ttide is a cor- 
rection factor that accounts for the fact that mass only needs 
to reach the Hill radius to escape (Erkaev et al. 2007). For 
planets like Kepler- 1 lb today this correction factor increases 
the mass loss rate by ~ 10%, however at early times it can 
increase the rate by as much as a factor of 2. 

Finally, e is an efficiency factor that parametrizes the frac- 
tion of the incident XUV flux that is converted into us- 
able work. This efficiency is set by radiative cooling, es- 
pecially via Ly man a, and can depe nd on the level of in- 
cident flux (M urrav-Clav et al.ll2009b . Kepler- 11 is a 8±2 
Gyr old Sun-like star. Using the power law from Ribas 
et al., this implies that current XUV flux at Kepler-llf 
is w 37 ergs"' cm"^. Similarly, when Kepler-11 was 100 
Myr old, the flux at Kepler-llb was w 6 x lO"' ergs"' cm"^. 
iMurrav-Clav et al.l (|2009) found that at XUV fluxes over 
10^ ergs"' cm"^, relevant for many hot Jupiters, mass loss be- 
comes radiation/recombination-limited and highly inefficient. 
However, at the lower XUV fluxes relevant for the Kepler- 1 1 
system mass loss is roughly linear with Fxuv and has efficien- 
cies ^ 0.1-0.3. For this work, we assume a default efficiency 
of e = 0.1±[|;i5, although we do examine the effects of lower 
efficiencies. While we predominantly investigate the loss of 
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H/He envelopes, in some limited cases for Kepler- lib, we 
also assume this holds for steam envelope loss. 

One important implication from equation (|2) is that mass 
loss rates are much higher when planets are young. This is 
due to two reasons. Planetary radii are considerably larger 
due to residual heat from formation. Moreover, at 100 Myr 
Fxuv was ~ 500 times higher than it is currently (IRibas et al.l 
l2005h . As a result, most of the mass loss happens in a planet's 
first Gyr. Thus although a planet's envelope may be stable to- 
day, its composition may have changed significantly since for- 
mation. Likewise, a considerable amount of mass wi ll be lost 
betwe en the end of planet formation at ^ 10 Myr ( Calvet et al.l 
20021) and 100 Myr . Following the x-ray observations of 



Jackson et al.l (1201 Ih . we assume that at ages younger than 



100 Myr the stellar XUV flux saturates and is constant at the 
100 Myr value. Unfortunately, the observations for 10-100 
Myr do not cover the EUV (100-1200 A) part of the spec- 
trum, so there is some uncertainty as to whether this saturation 
age is uniform across the entire XUV spectrum. Nonethel ess, 
our constraints on the formation of Kepler- 1 1 in sections 13.21 
and 15.11 come from the lower limits we are able to place on 
the initial compositions. Assuming that the EUV saturates 
along with the x-rays is conservative assumption in terms of 
the amount of mass that is lost. 

In general, models of LMLD planets that assume H/He en- 
velopes today will predict much larger mass loss histories then 
models that assume steam envelopes. Partly, this is because 
the lower mean molecular weight of hydrogen. Mostly, how- 
ever, it is because when we integrate the compositions back 
in time from the present, the addition of a small amount of 
H/He has much larger impact on a planet's radius than a small 
amount of water A larger radius in the past in turn means a 
higher mass loss rate; and so the integrated mass loss history 
becomes much more substantial for H/He envelopes. 

3 . APPLICATION TO KEPLER- 1 1 

3.L Current Compositions from Thermal Evolution 

The first step in trying to understand the formation and his- 
tory of a planetary system is to identify the possible current 
compositions for each of the planets in the absence of any 
mass loss. This then gives us estimates for the current masses 
of each planet's core, which we then use as the starting point 
for all of our calculations with mass loss. 

Figure [T] shows the Kepler- 11 planets in a mass-radius 
diagram along with curves for different possible com- 
positions. For all planets, we color-code by the inci- 
dent bolometric flux they receive. The Kepler- 11 plan- 
ets are shown by filled circles with identifying letters 
next to each one. The other known transiting exoplan- 
ets in this mass and radius range are shown by the open 
squares. In order o f increasing radi us, these are Kepler- 10b 
dBatalha etal."2011), Kepler-36b (Carter et al. 2012), CoRoT- 
7b (.Leger et al. 2009: Oueloz et al. 2009; Hat zes et al...201li 



Kepl er-20b dFressin et al.. 
18b dCochran etaLl 120111)" 



201lU Gautier et a l.ll2011h. Kepler- 
55 Cancri e (Winn et al.' '201 1^, 



iDemorvet all 1201 ik GJ 1214 b dCharb onneau et al. 2009), 
Kepler-36c (Carter etal. 2012'), Kepler-30b (Fabrvckv et alJ 
2012; Sanchis-Ojeda et al. 2012), and GJ 3470b (Bonfils et al] 
2OI2I) . Lastly, the open triangles show the four planets in our 
own solar system that fall in this range: Venus, Earth, Uranus, 
and Neptune. 

The curves show various possible compositions. The solid 
black curve shows a standard Earth-Hke composition with 2/3 



rock and 1/3 iron as described in Section l2~TI The other curves 
show compositions with thick water or H/He envelopes atop 
an Earth-like core. These curves include thermal evolution 
without mass loss to 8 Gyr, the age of Kepler-11. The blue 
dashed curves show the results for 50% and 100% water- 
worlds computed at Teq = 700 K, approximately the average 
temperature of the five inaner planets. Likewise, the dot- 
ted orange curves show the results for H/He envelopes; how- 
ever, here each curve is tailored to match a specific Kepler- 1 1 
planet and is computed at the flux of that planet. These fits 
are listed in greater detail in table [T] He r e we li st the mass 
of each planet taken from ILissauer et al.l (1201 lal) : the H/He 
fractions needed to match each planet's current radius for a 
water-poor super-Earth model; the water fraction needed to 
match Kepler- lib as a water- world; and the H/He fractions 
needed to fit Kepler- 1 Ic, d, e, and f as sub-Neptunes with an 
intermediate wat er lay er, as described in Section ED As de- 
scribed in section[32] we varied the planetary albedo, the heat 
capacity of the rocky core, and the observed mass, radius, cur- 
rent age, and incident flux. 

Figure [T] and table [T] clearly show the degeneracy between 
various compositions that we are attempting to untangle. 
There are now four planets including Kepler- 1 lb that can eas- 
ily be fit either as water-worlds or as water-rich sub-Neptunes 
with < 2% of their mass in H/He. However, it is worth look- 
ing closer at Kepler- 1 lb in particular It is this the only planet 
in the system which does not require any hydrogen or helium 
to match its current radius, although it must have some sort 
of volatile envelope. Moreover, it is also the most irradiated 
and it is fairly low gravity. As a result, adding a small amount 
of hydrogen to its current composition has a large impact on 
the bulk density, which in turn makes the planet more vul- 
nerable to mass loss, as seen in Eq. (|2]l. A clearer picture 
for this planet emerges when including XUV driven mass-loss 
and relatively str ong constraints from formation models dis- 
cussed in Section lsTl] Thus, if there is hope of using mass loss 
to constrain the composition and formation of the system, it 
likely lies with Kepler- 1 lb. 

3.2. Mass Loss for a Super-Earth Scenario 

Now that we have estimates for the present day compo- 
sitions, we will begin considering the effects of mass loss. 
We will compute mass loss histories that when evolved to the 
present day, match the current mass and composition. This 
then tells us what the mass would have to be in the past to 
result in the current mass and composition. As discussed in 
Section 12.31 there is uncertainty in stellar XUV fluxes ages 
younger than 100 Myr; as a result, we will present results 
both at 10 Myr and 100 Myr after planet formation. 

First we will consider water-poor super-Earth models for 
each planet, which have H/He envelopes atop Earth-like rocky 
cores. As discussed in Section l273l H/He envelopes are par- 
ticularly susceptible to mass loss. As an example. Figure |2] 
shows four possible cooling histories for Kepler- lib. The 
solid lines show thermal evolution without any mass loss 
while the dashed lines include mass loss. The orange curves 
are for water-poor super-Earth models, while the blue curves 
show water-world models. The red cross shows the current 
radius and age of Kepler- 1 lb. These curves illustrate the im- 
pacts of both thermal evolution and mass loss on the radius 
of a low-mass planet. The water-world models require that 
40% of the current mass must be in water to match the current 
radius. Assuming our standard efficiency e = 0.1, implies an 
initial composition of 43% water at 10 Myr. This illustrates 



5 



TABLE 1 
Current Mass and Composition 



Planet 


Current Mass (M^) 


% H/He 


% Water 


% H/He 3-layer 


Kepler-llb 




0.3±ii5% 


40 ±^1% 


n/a 


Kepler- 1 Ic 


13 5+''-* 


4.6±|]% 


n/a 


0.3±|!,'j'% 


Kepler-lld 


6.1±3.i 




n/a 


1.3±!]|% 


Kepler- 1 le 




17.2±^;l % 


n/a 




Kepler-llf 


2.3±ji 


4.1±|J% 


n/a 


0.4±|];«% 



Note. — Present day masses and compositions for Kepler-1 1 for three classes 
of models. The first two columns list each planet along with the observed mass 
derived from TTV in Lissauer et al. (2011a). The third column lists the cuiTent 
H/He fractions predicted by our theiTnal evolution models assuming a two layer 
H/He on rock/iron or water-poor "super-Earth" model. The fourth column shows 
the predicted water fraction for a two layer steam on rock/iron model or "water- 
world." This is only applicable to lib as the other are all less dense than pure 
water The final column lists the predicted fraction of H/He for three layer "sub- 
Neptune" models with equal mass in the rock and water layer This scenaiio is not 
applicable to 1 lb since it only needs 40% water to match its radius. 



TABLE 2 

Results from Mass Loss: Mass and Composition for Super-Earths 



Planet 


Mass 100 Myr (M®) 


% H/He 100 Myr 


Ma.ss 10 Myr (M^) 


% H/He 10 Myr 


Kepler-llb 




87.6 ±^514% 




90.4 ±|l% 


Kepler-1 Ic 


13.7±^;| 


6.0±5;0% 


14.2±j;3 


9.1 ±283% 


Kepler-lld 




16.5 ±^5% 




28 iff % 


Kepler-1 le 




21.2 ±5-°% 


9.7±2.5 


28.1 ±^"7^% 


Kepler-llf 


3.1±5;2 


29±f^% 


3A±%\ 


35 % 



Note. — Masses and H/He fractions predicted by coupled mass loss and thermal evolution 
models at 100 and 10 Myr, assuming all five planets are water-poor super-Earths. The large error 
bars on some compositions are due mostly to uncertainties in the current masses fr om T TV. The 10 
Myr values are subject to some model uncertainties as discussed in sections |T2l and |231 Kepler-1 lb 
is extremely vulnerable to H/He mass loss and would have to start off implausibly massive to retain 
a small H/He envelope today. For the best-fit masses, Kepler-1 Ic is less vulnerable to mass loss due 
to its massive core. However, Kepler-1 Ic, d, e, and f are all consistent with have foiTned with ~30% 
H/He. 



the relative stability of water envelopes. On the other hand, 
the dashed orange curve shows the vulnerability of H/He lay- 
ers. Here we have assumed a efficiency 5 x lower e = 0.02 and 
yet more mass is lost than in the water-world scenario. Even 
at this low efficiency, Kepler-llb would have to initially be 
1 1 % H/He and 4.8 to retain the 0.3% needed to match the 
current radius. This also shows the large increase in radius 
that can result from even a relatively modest increase in the 
H/He mass. 

Table |2] summarizes the results for Kepler- 11 b-f for the 
water-poor super-Earth scenario. We list the masses predicted 
by our models when the planets were 10 and 100 Myr old. 
In addition, we list the fraction of the planets' masses in the 
H/He envelope at each age. These results are further illus- 
trated in Figure H^. Here we have plotted the mass and H/He 
fraction for each planet at 10 Myr, 100 Myr, and today. Each 
color corresponds to a particular planet with the squares in- 
dicating the current masses and compositions, the circles the 



results at 100 Myr, and the triangles the results at 10 Myr. In 
order to calculate the uncertainty on these results, we varied 
the mass loss efficiency e from 5-20% and varied the XUV 
photosphere from 0.1-10 nbar Likewise, we varied the plan- 
etary Bond albedo from 0-0.80 and varied the heat capacity 
of the roc ky c ore from 0.5-1.0 Jg~'K~'. Also, as discussed 
in section 122] we varied the initial entropy for the 10 Myr 
models, to account for undertainties in "hot-start" vs. "cold- 
start". Finally, we factored in the observed uncertainties in 
mass, radius, and incident flux. 

Clearly, Kepler-llb is vulnerable to extreme mass loss if 
it has a H/He envelope atop a rock/iron core. Although less 
than 1 % H/He today, if it is a water-poor super-Earth it could 
have been have over ^ 90% H/He in the past. At 10 Myr, its 
mass would have been 45 ± 10 M®, an order of magnitude 
higher than the current value. Kepler-llb is able to undergo 
such extreme mass loss because its high XUV flux and the low 
mass of its rocky core put it in a regime where it is possible 
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to enter a type of runaway mass loss. This happens when the 
mass loss timescale is significantly shorter then the cooling 
timescale. After the planet initially loses mass it has an inte- 
rior adiabat and rocky core that are significantly hotter than 
would otherwise be expected for a planet of its mass and age. 
This is because the interior still remembers when the planet 
was more massive and has not had sufficient time to cool. As a 
result, the planet will stay inflated for some time and the den- 
sity stays roughly con stant and can actually decrease. A simi- 
lar effect was seen by Ba raffe et al.l (12004 ) when they studied 
coupled thermal evolution and mass loss models for core-less 
hot Jupiters. We find that this process generally shuts off once 
the composition drops below ^20% H/He. At that point the 
presence of the core forces the total radius to shrink even if 
the planet is unable to cool efficiently. Figure |3] shows this 
process as Kepler- 1 lb loses mass for three different values of 
its current mass and therefore its core mass. The curves cor- 
respond to the best fit mass from transit-timing as well as the 
la error bars. This shows that the timing of this runaway loss 
event depends strongly on the mass of the rock/iron core. 

Super-Earth models of Kepler- 1 lb are unusual in that they 
are subject to tremendous mass loss and yet they retain a small 
amount of H/He today. Typically models that start out ^ 90% 
H/He either experience runaway mass loss and lose their H/He 
envelopes completely, or they never enter the runaway regime 
and remain over 50% H/He. The uncertainty in the initial 
composition of Kepler- lib is due to uncertainty in its TTV 
mass. At a given current mass, the range of Kepler- lib mod- 
els that will retain an envelope that is < 1 % H/He is extremely 
narrow. In this sense, the current composition of Kepler-1 lb 
requires a rare set of initial conditions if it is a water-poor 
super-Earth. 

Counterintuitively, if Kepler-1 lb is more massive today 
then its implied mass in the past is actually lower. This is be- 
cause a higher mass today would imply a more massive core, 
which would increase the planet's density and decrease its 
mass loss rate. As a result, a more massive model for Kepler- 
1 lb today is less vulnerable to mass loss and so less H/He is 
needed in the past in order to retain 0.3% today. At 100 Myr, 
there is a very large uncertainty in the composition due to the 
uncertainty in the core mass. However, even if we assume the 
Icr error bar 6.5 M^, Kepler-1 lb would still be at least 37 
and at least 83% H/He at 10 Myr In section lSTTl we will com- 
pare this to models of in situ formation and show that such a 
scenario is unlikely. 

On the other hand, Kepler- 1 1 c is not particularly vulnerable 
to mass loss, at least using the best fit mass from transit tim- 
ing, despite having the second highest flux in system. This is 
because of the relatively large mass of its rocky core; the high 
gravity means additional H/He has a more modest effect of 
the planet's radius and therefore on the mass loss rate. In fact, 
along with the incident XUV flux the mass of the rocky core 
is the single largest factor that determines whether a given 
planet will be vulnerable to mass loss. As a result, the dom- 
inant sources of uncertainty in our mass loss models are the 
uncertainties in the masses from TTV. These dominate over 
all the theoretical uncertainties in the thermal evolution and 
mass loss models. The uncertainty in planet mass from transit 
timing is particular large for Kepler-1 Ic. If its mass is close 
to the la low value, then it is possible Kepler-1 Ic has un- 
dergone more substantial mass loss similar to Kepler-1 Id-f. 
Fortunately, as more quarters of data are processe d the mass 
estimates from TTV will become more precise ("Ag ol et all 
12005; .Holman & Murray .2005.) . Finally, Kepler-1 Id, e, and 



f are modestly vulnerable to mass loss and are consistent with 
having originated with 20% H/He at 100 Myr and ~ 30% 
H/He at 10 Myr In Section l572l we will discuss these results 
in terms of orbital stability. 

3.3. The Water-Rich Scenario 

Next we consider a water-rich scenario where the entire 
system formed beyond the snow line. We assume that Ke pler- 
1 Ic-f are water-rich sub-Neptunes as described in Section ETTl 
while Kepler-1 lb is currently a water- world. Otherwise the 
thermal mass loss histories are calculated in the same man- 
ner as the water-poor super-Earth scenario. For Kepler-1 Ic- 
f we calculate the planet mass H/He fraction at 10 and 100 
Myr, assuming that only H/He is lost. For Kepler-1 lb, we 
examine the vulnerability of both H/He and steam envelopes 
atop water-rich interiors. The results are summarized in Ta- 
ble[3]which list the water fraction for a water-world model of 
Kepler-1 lb and the H/He fraction for water-rich sub-Neptune 
models of Kepler-1 Ic-f. Likewise, the results for c-f are 
shown in Figure l^^- 

In general, these three layer models are slightly more vul- 
nerable to mass loss than the water-poor super-Earth models 
presented in section [J!2] Mostly this is because models with a 
water layer have hotter interiors that cool more slowly. Since 
models For example, for Kepler-1 Ic without mass loss the 
models presented in Table [T] at 8 Gyr the final entropy in the 
H/He layer is 6.6 kb/baryon for the water-rich sub-Neptune 
model versus 5.8 for the water-poor super-Earth model. The 
second reason is that counter-intuitively the water-rich sub- 
Neptune models are slightly more vulnerable to mass loss pre- 
cisely because they have less of the planet's mass in H/He to- 
day. For a planet that has less H/He today, adding a small 
amount of H/He at the margin has a larger impact on the 
planet's radius and therefore on the mass loss rate. 

For Kepler-1 Ic-f the results are broadly similar to the those 
for the water-poor super-Earth scenario. Kepler-1 Ic is again 
the least vulnerable to mass loss; while Kepler-1 Id is again 
the most vulnerable of the four planets that we model as 
water-rich sub-Neptunes. However, all four of these planets 
are consistent with having been ^ 10-20% H/He at 100 Myr 
and - 20-30% H/He at 10 Myr 

If Kepler-1 lb was always a water- world, then mass loss 
was never important for it. Between 10 Myr and the present, 
it only drops from 43% to 40% water Moreover, if Kepler- 
1 lb was initially a water-rich sub-Neptune similar to the other 
planets in the system, it could have easily stripped its H/He 
outer envelope. If we start Kepler-1 lb at 100 Myr as a water- 
rich sub-Neptune similar to the other planets with 30% H/He 
atop 4.3 M0 of rock and water, then assuming e = 0.1 the 
entire H/He envelope will be stripped by 300 Myr We can 
set upper limits on the initial mass and H/He fraction of 70 
M© and 94% if Kepler-1 lb was originally a water-rich sub- 
Neptune. These are however strictly upper limits, a H/He 
layer could have been lost at any time between formation and 
now. Therefore, all five planets are consistent with a scenario 
is which they formed as water-rich sub-Neptunes with ^ 10% 
H/He at 100 Myr and - 20% H/He at 10 Myr 

4. A MASS LOSS THRESHOLD FOR LOW-MASS LOW-DENSITY 
PLANETS 

Although Kepler- 11 provides a unique case-study, it is es- 
sential to explore how mass loss impacts the larger population 
of LMLD transiting planets. Figure |5] shows the bolomet- 
ric flux these planets receive at the top of their atmospheres 
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TABLE 3 

Results from Mass Loss: Water-Worlds and Sub-Neptunes 



Planet Mass 100 Myr (Mq) Composition 100 Myr Mass 10 Myr (M®) Composition 10 Myr 



Kepler- lib 




41 ±11 % 


4 5±2' 


43 ±38 % 


Kepler- 1 Ic 


13 7+"-* 


1.8 ±1*4% 


15.2±f2 


12±™% 


Kepler- lid 


6.8±2;8 


1 1.6 ±1^7% 


7.6±^;« 


21±|9% 


Kepler- 1 le 




12.8 ±^-5% 


9.7±2;« 


18±|2% 


Kepler-llf 




21 ±« % 


4.0±|;« 


43 ±11 % 



Note. — Masses and volatile fractions predicted by coupled mass loss and thermal evolution models at 
100 and 10 Myr. assuming Kepler-llc, d, e, and f are water-rich sub-Neptunes and Kepler-llb is a water- 
world. Thus the compositions listed for Kepler- 1 lb are water fractions, while those for Kepler- 1 Ic-f are H/He 
fractions. Kepler- 1 Ic-f are all consistent with having formed as water-rich sub-Neptunes with 20-30% H/He. 
Kepler- 1 lb is not vulnerable to mass loss if it has a water envelope; however, it could have easily also formed 
a water-rich sub-Neptune and lost its H/He envelope. 



VS. their bulk densities. As in Figure [T] filled circles show 
the Kepler- 11 planets with the letters indicating each planet. 
Likewise, the open squares show the other transiting exoplan- 
ets that are less than 15 Earth masses. For reference, we have 
also plotted all other transiting planets between 15 and 100 
M0 as gray crosses (Wright et al. 2011). The colors indi- 
cate possible compositions. All planets with a best-fit mass 
and radius that lies below a pure rock curve are colored red. 
These include Kepler-lOb, Kepler-36b, CoRoT-7b, and just 
barely Kepler-20b. Planets that are less dense than pure rock 
but more dense than pure water, indicating that the could po- 
tentially be water-worlds, are colored blue. These include 
Kepler-llb, Kepler-18b, and 55 Cancri e. Meanwhile those 
planets that must have a H/He envelope to match their radius 
are colored orange. These include Kepler-llc, d, e, and f, 
Kepler-30b, Kepler-36c, GJ 1214b, and GJ 3470b. 

The dashed black lines show curves of constant mass loss 
rate according to equation assuming e = 0.1 and /iTtide = 1- 
These curves are linear in this plot since the instantaneous 
mass loss rate goes as the flux over the density. Although Fig- 
ure |5] plots the bolometric flux toda y, we can re l ate th is to an 
XUV flux at a given time using the iRibas et aTl (l2005h power 
law for sun-like stars described in Section 12.31 The curves 
show the flux today required to lose mass at 1 Gyr"' when 
the planets were 1 Gyr old and 100 Myr old, along with an- 
other curve showing 0.1 Gyr"' at 100 Myr Since most of 
the mass loss happens in the first few hundred Myrs, the bot- 
tom two curves can roughly be considered as the respective 
thresholds for mass loss being important and being unimpor- 
tant for LMLD planets. 

One possible explanation of this mass loss threshold is that 
it caused by XUV driven mass loss from H/He envelopes on 
low-mass planets. LMLD planets that form above the 100 
Myr 1 M^Gyr"' curve lose mass, increase in density and 
move to the right until they lie below this threshold. The 
planets that are left above this line are mostly rocky or at 
the very least probably do not have H/He envelopes. Planets 
more massive than ~ 15 are not affected since they have 
a larger reservoir of mass and the loss of a few earth masses 
of volatiles isn't sufficient to significantly change their bulk 
density. To illustrate this, we have plotted our predictions for 
the bulk densities of each of the Kepler- 1 1 planets at 100 Myr, 
including the effects of both mass loss and thermal evolution. 



These are indicated by the shadowed letters at the left of Fig- 
ure|5] 

The situation becomes even clearer if we instead we plot 
flux against mass times density as in Figure |6] The timescale 
for XUV mass loss goes like pMp/Fxvv, so lines in this di- 
agram are constant mass loss timescales. Now the thresh- 
old is much clearer and applies to all planets up to all plan- 
ets with H/He envelopes. This also removes any effects from 
the somewhat arbitrary 15 cut. The sparsity of planets at 
low flux and high density is almost certainly a selection effect, 
since these are likely to be planets with long periods and small 
radii. However, the interesting result is that there is appears to 
be a critical mass loss timescale above which we do not find 
any planets with H/He envelopes. In particular, all five of the 
inner Kepler- 1 1 planets lie nicely along this threshold. More- 
over, of the three planets that lie above the critical mass loss 
timescale, two are likely rocky. 



M 7re/?^Fxuv.Ei()0 F^, 

The dashed black line in Figure |6] shows our best fit for 
this critical mass loss timescale. Equation |5] defines this 
mass loss timescale. Here e = 0.1 is the mass loss efficiency, 
■Pxuv.Eioo = 504 erg s"' cm"^ is the XUV flux at the Earth when 
it was 100 Myr old, and Fp is the current incident bolometric 
flux at a planet. We find a best fit with rio,ss,crit ~ 12 Gyr How- 
ever, while equation |5] accounts for the higher XUV fluxes at 
earlier times, it does not include the effects of larger radii at 
formation. The will reduce fioss by at least another order of 
magnitude. 

A similar mass l oss th reshold was proposed by 
Lecavelier Pes Etangsl (I2007h . Unfortunately, at that 
time there were relatively few transiting planets and no 
known transiting super-Earths. As a result, the authors we 
mostly limited to hot Jupiters from radial velocity surveys 
and were forced to use a scaling law to estimate radii. Here 
we are able to confirm the existence of a mass loss threshold 
and extend it all the way down to ^ 2 M^. 

This mass loss threshold could also help ex plain features 
in occ urrence rate of planets found by Kepler. [Howard et al.l 
( I2OI Ih found that the frequency of 2-4 Kepler planet can- 
didates dropped off exponentially for periods within 7 days. 
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This 7 day cutoff corresponds to an incident bolometric flux 
of 200 Fq. There are five planets with measured densities in 
figure |5] that lie above 200 F^. Of these five, three planets 
are consistent with being rocky and two with being water- 
worlds; none of the five requires a H/He atmosphere to match 
its observed mass and radius. If all low mass planets orbiting 
within 7 days lose their H/He atmospheres, then their radii 
will shrink from 2-4 to <2 Rq. This could naturally ex- 
plain the drop off in 2-4 R^ candidates at short periods. 



4.1. Reproducing the Mass Loss Threshold 

In order to fully examine whether the mass loss threshold in 
Figure |5] can be explained by atmospheric mass loss, we per- 
formed a small parameter study with ^800 mass loss models 
across a wide range of initial masses, compositions, and in- 
cident fluxes. For each model we ran thermal evolution and 
mass loss starting at 10 Myr around a Sun-like star We ran 
models with initial masses of 2, 4, 8, 16, 32, and 64 Mq. 
We assumed water-poor super-Earth compositions, meaning 
H/He envelopes on top Earth-like cores, with initial compo- 
sitions of 1, 2, 5, 10, 20, and 40% H/He. Finally we varied 
the incident bolometric flux from 10 to 1000 Fm, in order to 
cover the range of observed planets in Figures js] and |6] We 
then recorded the resulting masses, densities, and composi- 
tions at various ages. 

The results are shown in Figure |7] As in Figure |6] each 
panel plots the total incident flux at the top of the atmosphere 
vs. the planet mass times density assuming different mass loss 
histories for our full suite of models. The size of each point 
indicates the mass of the planet, while the color indicates the 
fraction of its mass in the H/He envelope. The top left panel 
shows the initial distribution at 10 Myr before we start any 
mass loss. The other two top panels show the results at 100 
Myr and 10 Gyr for our standard mass loss efficiency e = 0.1. 
Meanwhile, the bottom panels show the results at 1 Gyr for a 
range of different efficiencies. These range from highly inef- 
ficient mass loss e = 0.01, to our standard efficiency e = 0.1, 
and finally extremely efficient mass loss e = 1 . In each panel, 
as planets cool and lose mass the points move to the right, 
shrink, and become bluer (less H/He). For reference, we have 
re-plotted our critical mass loss timescale from Figure 0in 
each of the result panels. 

As we can see, models with mass loss do in general re- 
sult in a threshold roughly corresponding to a critical mass 
loss timescale. Moreover, the mass loss threshold observed in 
Figure|5]is well reproduced by mass loss models with e w 0.1. 
This is similar to the efficiencies found by detailed models 
of mass loss from hot Ju piters in the energy-limited regime 
dMurrav-Clav et aLll2009h . This suggests that our assumption 
of comparable mass loss efficiencies for LMLD planets is rea- 
sonable. It is also apparent that the threshold already in place 
by 100 Myr, and subsequent evolution has a relatively minor 
effect. We also examined the effect of beginning our param- 
eter study at 100 rather than 10 Myr; however, this did not 
significantly affect the location of the threshold; 

Previous mass loss evolution models (e.g.. lHubbardetal.1 
I2007blat iJackson et al.) l20TTt lOwen & JacksonI l2012h have 
also predicted mass loss thresholds. However, our models are 
the first to fully include the effects of coupled mass loss and 
thermal evolution for LMLD planets. We are able confirm and 
explain the observed threshold seen in Figure |6]in a region of 
parameter space where most of the Kepler planets are being 
found. 



4.2. Constraints On Mass and Radius for the General 
Population 

If we use the critical mass loss timescale curve from Figure 
|6] as a approximation for the observed mass loss threshold, 
then we can write down a simple expression for the thresh- 
old. This is shown in equation which is valid for plan- 
ets around Sun-like stars with Fp < 500 F®. The 500 cut 
excludes highly irradiated rocky planets like Kepler- 10b and 
CoRoT-7b. These planets may have once had volatiles in the 
past, but they are likely rocky today and so H/He mass loss 
is no longer relevant. This cut also excludes the region where 
energy-limited escape breaks down and mass loss becomes ra - 
diation and recombination limited (IMurray-Clay et afl2009l) 



pMp> 



3eF. 



XUV.EIOO 



4G 



(6) 



The exciting implication of equation ^ is that we can use it 
to obtain lower limits on mass for the much larger population 
of Kepler super-Earths and sub-Neptunes for which we do not 
have measured densities. This will help identify promising 
targets for follow-up work with radial velocity observations. 
This is shown in equation (|7|. 
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Table |4] applies equation (|7]) to a list of Kepler candi- 
dates smaller than 4 R^ that are well suited to radial-velocity 
follow-up. We excluded any planets with fp > 500 F^, since 
equation (|7]) is not valid in that regime. Also, we limited 
the sample to only those planets with minimum radial veloc- 
ity semi-amplitudes K^nin > 1.0 ms"' around stars with Ke- 
pler magnitude brighter than 13, since these will be the most 
promising for RV follow-up. In the end, this leaves us with 
a list of 38 likely detectable targets, eight of which (KOIs 
104.01, 107.01, 123.01, 246.01, 262.02, 288.01, 984.01, and 
1241.02) have /Tmin > 2.0ms-i. 

Finally, we can also use the mass loss threshold to find an 
upper limit on the radii of non-transiting planets from radial 
velocity surveys with Fp < 500 F^. This is done in equation 
©. 
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5. DISCUSSION 



5.1. Kepler-11: Comparison to Formation Models, 
Implications for Migration 

By itself, the constraints from mass loss do not tell us 
whether Kepler-11 is a system of water-poor super-Earths or 
water-rich sub-Neptunes. Instead we need to compare our es- 
timates of the initial compositions to models of planet forma- 
tion. By doing so we can examine whether our estimates of 
the original compositions for a water-poor super-Earth sce- 
nario are consistent with the maximum H/He fraction that can 
be accreted during in situ formation. 

Ikoma & Hori (2012) examine the accretion of H/He atmo- 
spheres onto the rocky cores of hot water-poor super-Earths . 
In particular, they examine the in situ formation of the Kepler- 
1 1 system. In addition to a planet's core mass and tempera- 
ture, the amount of H/He accreted will depend strongly on the 
lifetime and dust grain opacity of the accretion disk. As with 
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thermal evolution, the need to cool the rocky core can slow 
the contraction of the accreting atmosphere and limit the final 
H/He fraction. They are able to set hard upper limits on the 
initial compositions for in situ formation by assuming a grain- 
free, long-lived 1 Myr in the inner 0.2 AU) accretion disk 
and ignoring the delay in accret i on due to cooling the core. 

In paiticulai-, llkoma & Horil (120121) find that Kepler- lib 
could not have accreted more than 10% of its mass in H/He 
if it formed in situ. Moreover, using a more typical disk life- 
time of 10^ yr (t Gorti et all 120091) and including the effect of 
cooling the core implies that Kepler- lib was < 1% H/He at 
formation. On the other hand, in Table|2]we showed that ther- 
mal evolution and mass loss models predict that if Kepler- 1 lb 
is a water-poor super-Earth then it was 87 and at least 

82% at 10 Myr Combined with the results of 'Ikoma & Horil 
(j2012), this disfavors in situ formation of Kepler-llb. This 
result appears robust to any uncertainties in thermal evolution 
or mass loss models. Even if we only look after the period of 
run-away mass loss, at 3 Gyr Kepler- 1 lb was s till 10% H/He, 
the maximum allowed bv llkoma & Horil (120 121) . Likewise, we 
find that K epler- 1 If was at leas t 10% H/He at 10 Myr, even 
though the Ikoma & Horil ('2012) models predict that it cannot 
have accreted its current composition of 4% H/He if it formed 
in situ. Furthermore, the co-planar, tightly packed, circular 
orbits in the system stron gly suggest that it could have un- 
dergone type 1 migration (llda & LinI 120101) . As a result, we 
disfavor in situ formation of the system. 

If the Kepler- 11 system did not form at its current loca- 
tion, then one possibility is that it formed at or beyond the 
snow-line and then Type 1 migrated to it is current location 
dRogers et al.l 1201 H) . If this is the case, then it is likely a 
system of water-rich sub-Neptunes and water-w orld s as dis- 
cussed in Section [33] As we showed in Section [331 Kepler- 
1 lb is very stable to mass-loss if it is a water-world. If it was 
initially a water-rich sub-Neptune, it could have easily lost its 
H/He layer in the first few 100 Myr Likewise, Kepler- llc-f 
are all consistent having formed as water-rich sub-Neptunes 
with ~ 20% of their mass in H/He. 

The other possibility is that Kepler- 1 1 is a system of water- 
poor super-Earths that has nonetheless undergone significant 
migration. For a grain-free accretion disk that lasts 10^ yr 
at 5 50 K, the critical mas s for run-away accretion drops to 5 
(llkoma & Horil2012l) . This impHes that Kepler- 1 lb could 
possibly have formed as a water-poor super-Earth at or be- 
yond the current orbit of llf. Nonetheless, this assumes a 
completely grain-free long-lived disk, which may not be real- 
istic. Furthermore, this scenario still requires that Kepler-1 lb 
was ^90% H/He when it formed, while all the other planets 
in the system are consistent with more modest initial com- 
positions. As a result, we favor the water-rich sub-Neptune 
scenario. 

5.2. Kepler 11: Mass Loss and Orbital Stability 

One possible result of significant mass loss is that it could 
impact the orbital stability of closely packed multi-planet sys- 
tems like Kepler-1 1. Although this system is stable in its cur- 
rent configuration, it might not be with the initial masses de- 
termined by our models. One relatively simple stability check 
is to calculate the separation between pairs of planets in terms 
of their mutual Hill spheres. Figure |8] plots the separation in 
mutual Hill spheres (A) between adjacent pairs of planets at 
both 10 Myr and the present, assuming a water-poor super- 
Earth composition. 



ISmith & Lissauerl (120091) found that systems with five or 
more planets tended to de-stabilize when A < 9. This 
threshold is shown as dashed gray lines in figure |8] Al- 
though Kepler - 1 Ib-c c urrently lies well below this threshold, 
Lissau er et al.l (1201 lal) showed that the system is nevertheless 
stable today because planets b and c are dynamically decou- 
pled from the other four planets and so act more like a two 
planet system. For two planet systems the absolute minimum 
stable separation is A = 2^3 = 3.46 ( lGladmanlll993l) . This 
second stability threshold is shown by the dotted lines in fig- 

ure[8] The Hill radius goes as Mp^^, as a result the change in A 
from mass loss is relatively modest; nonetheless, the stability 
of the system is in danger At 10 Myr, Planets d-e do lie be- 
low the approximate A > 9 stability threshold; however, both 
pairs on either side of d-e are still relatively stable which may 
help stabilize the system. More importantly, the separation of 
planets b-c at 10 Myr skirts dangerously close, A = 3.8±q 4, to 
the critical A > 2\/3 stability threshold. More detailed mod- 
eling needs to be done to assess the impact of mass loss on 
orbital stability; nonetheless, the A > 2\/3 stability threshold 
provides another strong reason to be skeptical of a water-poor 
super-Earth scenario for Kepler-1 lb. 

The major caveat to this stability analysis is that we assume 
that all of the orbits are s tationary even as the planets lose 
mass. This is motivated by lAdamsr(l201 ll) . which showed that 
in the presence of a modest planetary magnetic field XUV 
driven mass loss from hot Jupiters tends to come out along 
the magnetic poles. Assuming that the magnetic field is suf- 
ficiently strong, dipolar, and perpendicular to the plane of the 
orbit, then mass loss won't have any impact on the orbit. In 
general however, the directionality of mass loss will be an ex- 
tremely complicated problem determined by the interaction of 
the ionized hydrodynamic wind, the planetary magnetic field, 
and the stellar wind. Boue et al. (2012) showed that if the 
mass loss is directed in the plane of the orbit, then it can have 
a significant impact on both semi-major axis and eccentricity. 

5.3. Future Work 

Further confirmation of the mass loss threshold will de- 
pend on getting reliable mass estimates for more LMLD 
planets. Fortunately, there is a large population of super- 
Earth sized planets in multi-planet systems found by Kepler 
dLissauer et al.ll201 IH) . For so me of these syste ms TTV can 
be used to determine masses dAgol et all 120051) . Moreover, 
this will become possible for more systems as more quarters 
of data are collected. Likewise, as more quarters of transit 
data are analyzed, previous mass constraints from TTV will 
become more precise. For Kepler- 11 this will allow tighter 
constraints on both the current and past compositions. 

In order to better understand mass loss, there is also a 
strong need to acquire more XUV observations of young 
10 Myr) stars. Currently the best estimat es of EUV 
fluxes are for G and K stars o lder than 100 (iRibas et al.l 
120051; ISanz-Forcada et all I2010I) . Meanwhile, planet for- 
matio n ends and mass l o ss becomes important after a few 
Myr dCalvet et al.l 120021; [Alexander et"al] 120061) . A large 
amount of planetary mass will be lost in the first 100 Myr 
and this will depen d strongly on the stellar XUV flux. 
iJackson et al.l d201 11) recently found that x-ray fluxes satu- 
rate for stars younger than 100 Myr; however, this needs 
to be investigated at other wavelengths. Likewise, more 
UV observations of transiting exospheres are needed. Cur- 
rently, we only have observations for a handful of hot Jupiters 
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dVidal-M adiar et al.|l2004l: iLecaveUer Pes Etangs et al.ll2010l: 
iLecavel ier des Etangs et alj |2012|k More observations are 
needed, especially for Neptune and super-Earth sized planets. 

Likewise, it is also important get more observations of 
XUV fluxes and flares from M dwarfs, which are known to be 
highly active (Reiners et al. 2012). This is particularly impor- 
tant for super-Earths in the habitable zones of late M dwarfs 
where the stability of habitable atmospheres could depend on 
XUV driven mass loss. The three planets orbiting M-stars 
in figure |6] all lie an order of magnitude below the threshold 
for Sun-like stars. This could be due to mass loss from X- 
ray flares, but wit hout a larger s ample size it is impossible to 
say. The Recentlv lFrance et al.l ( 120121) obtained the first FUV 
spectrum for GJ876; however, more data are needed. Like- 
wise, non-equilibrium mass loss models need to be developed 
to understand the impact of flares. 

In the area of modeling mass loss, we have uti- 
lized a relatively simple model that appears to work 
well for hot Jupiters, although mas s loss rates have 
only been constrained for HD 209458b ("Vidal-Madi ar et al.l 
2004 and HP 18973 3b (LecaveHer Pes Etangs et al.i 120101: 
Lecavelier des Etangs et al.ll20 12). On a broader scale, there 
may be evidence for a lost population of h ot Jupiter planets at 
very h igh levels of XUV irradiation (e.g., iPavis & WheatlevI 
120091) . Since we have suggested that H/He mass loss may be 
more important for more modestly irradiated LMLP planets 
than for average hot Jupiters (since the smaller planets have 
much smaller H/He masses, mass loss for LMLP planets can 
actually transform the very structure of the planets), we en- 
courage further detailed models to quantitatively access mass 
loss from these atmospheres. 

6. CONCLUSIONS 

In order to better understand the structure, history, and for- 
mation of low-mass planets, we constructed coupled thermal 
evolution and mass loss models of water-poor super-Earths, 
water- worlds, and water-rich sub-Neptunes. The Kepler- 11 
system represents a new class of low-mass low-density plan- 
ets that offers a unique test-bed for such models and gives 
us powerful insights on planet formation and evolution. Ap- 
plying this understanding more broadly, we find a relation 
between a planet's mass, density, and its incident flux that 
matches the observed population. Moreover, this threshold 
can help constrain the properties of hundreds of planets. Our 
primary conclusions are: 



• XUV-driven hydrogen mass loss coupled with plane- 
tary thermal evolution is a powerful tool in understand- 
ing the composition and formation of low-mass low- 
density planets. 

• A coupled model is essential for this work, due to the 
much larger planetary radii in the past, when XUV 
fluxes were significantly higher 

• In situ formation of the Kepler- 1 1 system is disfavored, 
instead it could be a system of water-rich sub-Neptunes 
that formed beyond the snow line. 

• If Kepler- 1 lb is a water-poor super-Earth then it likely 
formed with - 90% H/He beyond 0.25 AU. We beUeve 
this is unlikely and instead show that Kepler- 11 b-f all 
could have originated as water-rich sub-Neptunes with 

20% H/He initially. If this is the case, Kepler- lib 
could have lost its H/He envelope and become a water- 
world today for a wide range of initial masses and com- 
positions. 

• There is a sharp observed threshold in incident flux 
vs. planet density times mass above which we do not 
find planets with H/He envelopes. To date, low-density 
planets have not been found above this threshold. 

• This mass loss threshold is well reproduced by our cou- 
pled thermal evolution and mass loss models. 

• This threshold can be used to provide limits on planet 
mass or radius for the large population of low-mass 
low-density planets without measured densities. 

• In particular, we have identified promising Kepler tar- 
gets for RV follow-up. 
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TABLE 4 

Minimum Masses and Densities for Selected KOI Candidates 



KOI 


KepMag 


P (day) 


Fp (Fe) 




Pmin (gem 3) 




Kmin (m s ') 


70.01 


12.50 


10.85 


80.3 


3.09 


0.73 


3.9 


1.21 


70.02 


12.50 


3.69 


343.5 


1.92 


3.0 


3.9 


1.75 


85.01 


11.02 


5.85 


403.4 


2.35 


2.4 


5.8 


1.83 


94.02 


12.21 


10.42 


209.9 


3.43 


1.0 


7.4 


1.92 


104.01 


12.90 


2.50 


233.4 


3.36 


1.1 


7.5 


4.08 


105.01 


12.87 


8.98 


130.3 


3.35 


0.82 


5.6 


1.94 


107.01 


12.70 


7.25 


301.7 


3.09 


1.4 


7.6 


2.29 


110.01 


12.66 


9.94 


220.5 


2.92 


1.3 


5.9 


1.62 


115.02 


12.79 


7.12 


409.1 


1.88 


3.4 


4.2 


1.29 


117.02 


12.49 


4.90 


436.9 


1.70 


4.1 


3.7 


1.25 


122.01 


12.35 


11.52 


108.8 


2.78 


1.0 


3.9 


1.06 


123.01 


12.37 


6.48 


461.4 


2.64 


2.2 


7.4 


2.46 


124.01 


12.94 


12.69 


227.8 


3.00 


1.2 


6.3 


1.65 


246.01 


10.00 


5.39 


404.8 


2.53 


2.2 


6.5 


2.27 


257.01 


10.87 


6.88 


308.6 


2.61 


1.8 


5.9 


1.80 


262.02 


10.42 


9.37 


491.8 


2.79 


2.1 


8.3 


2.20 


277.01 


11.87 


16.23 


177.4 


3.82 


0.79 


8.0 


1.91 


280.01 


11.07 


11.87 


154.9 


2.52 


1.3 


4.0 


1.16 


281.01 


11.95 


19.55 


192.3 


3.46 


0.95 


7.2 


1.99 
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TABLE 4 — Continued 



KOI 


KepMag 


P (day) 


(Fe) 




Pniin (g cm ^ ) 


Mp.min (Me) 


^min (ms ') 


285.01 


11.57 


13.74 


180.5 


3.38 


0.96 


6.7 


1.61 


288.01 


1 1.02 


10.27 


433.9 


3.11 


1.6 


9.2 


2.14 


291.02 


12.85 


8.12 


247.8 


2.14 


2.2 


3.9 


1.25 


295.01 


12.32 


5.31 


339.7 


1.77 


3.4 


3.5 


1.22 


297.01 


12.18 


5.65 


482.1 


1.65 


4.6 


3.7 


1.25 


301.01 


12.73 


6.00 


399.2 


1.75 


3.8 


3.7 


1.17 


323.01 


12.47 


5.83 


166.2 


2.17 


1.7 


3.3 


1.21 


984.01 


1 1.63 


4.28 


259.7 


3.19 


1.2 


7.4 


2.91 


987.01 


12.55 


3.17 


404.8 


1.28 


6.1 


2.3 


1.03 


1117.01 


12.81 


11.08 


327.5 


2.20 


2.4 


4.7 


1.13 


1220.01 


12.99 


6.40 


441.4 


1.95 


3.4 


4.6 


1.52 


1241.02 


12.44 


10.50 


485.3 


3.84 


1.3 


13.3 


3.17 


1597.01 


12.68 


7.79 


423.6 


2.67 


2.0 


7.2 


1.86 


1692.01 


12.56 


5.96 


175.9 


2.65 


1.3 


4.6 


1.61 


1781.01 


12.23 


7.83 


61.3 


3.29 


0.58 


3.7 


1.38 


1781.02 


12.23 


3.00 


219.6 


1.94 


2.4 


3.2 


1.63 


1921.01 


12.82 


16.00 


172.9 


3.09 


1.0 


5.7 


1.28 


1929.01 


12.73 


9.69 


251.7 


2.00 


2.4 


3.6 


1.11 


2067.01 


12.58 


13.24 


347.2 


2.97 


1.6 


7.6 


1.69 



Note. — Minimum masses and densities given by equation Q for KOI candidates selected for being promising targets for RV follow-up. In order to focus 
on promising candidates we cut the sample to planets with K^^in > 1.0ms"' around stars with KepMag < 13. We only included planets with Fp < 500 and 
Sp < 4 Sg). This leaves us with 38 high-quality targets, eight of which have K^in > 2.0 ms"'. 
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Planet Mass (Me) 

Fig. 1 . — Radius vs. mass for transiting exoplanets with measured masses, along witli curves for different compositions. Planets are color-coded by the incident 
bolometric flux they receive. Kepler-11 planets are shown by filled circles with letters indicating each planet. Other known exoplanets in this mass and radius 
range are shown by open squares. Solar system planets Earth, Venus. Uranus, and Neptune are shown by black letters. The solid black curve is for a Earth-like 
composition with 2/3 rock and 1/3 iron. All other curves use full thermal evolution calculations, assuming a volatile envelope atop a earth like core. The dashed 
blue curves are for 50% and 100% water by mass. The dotted orange curves are for H/He eiwelopes at 8 Gyr; each one is tailored to match a Kepler-1 1 planet 
and is computed at the appropiiate flux and for that planet. 
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Fig. 2. — Radius vs. time for four example model runs that match the present day mass and radius of Kepler-1 lb. The blue curves show water-world models, 
while the orange curves show water-poor super-Earth models. Dashed lines are with mass loss, while solid are without. Both water-world models and the 
water-poor super-Earth model without mass loss show very similar cooling curves. Even with our standard efficiency of 10% the water models undergo only 
minor mass loss. Meanwhile, even with an efficiency 4x smaller the H/He model undergoes substantial mass loss. This model is initially 14% H/He and 5.0 
Mfg. We have marked the masses for the H/He with mass loss model at 10 Myr, 100 Myr, and 1 Gyr This also shows the large impact on radius that even a 
modest (compared to Figure|3] H/He envelope can have. 
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Fig. 3 . — Mass vs. time with mass loss for three model mns that match the present day mass along with its Icr range for Kepler-1 lb. All three models assume 
a water-poor super-Earth composition that is 0.3% H/He today. The curves are color-coded by log density. The solid line corresponds to the best fit current mass 
from TTV; the dashed lines correspond to the 1 a bounds. This demonstrates several features described in the text. The initial mass is actually lower if Kepler- lib 
is more massive today due to a correspondingly more massive core. There is a period of runaway mass loss duiing which the density actually declines slightly, 
and the timing of this period depends strongly on the mass of the rocky core. 
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Fig. 4. — Composition vs. mass for models of Kepler-1 1 witli mass loss. Panel a) shows the results for water-poor super-Earth models with a H/He envelope 
atop a rock/iron core. Panel b) shows the results for water-rich sub-Neptune models that also have a thick water layer in between. Each point shows the % 
H/He and mass predicted by our thermal evolution and mass loss models at a given time. Each color indicates a particular planet as identified by the letters and 
connected by dashed lines. The open squares show the present day mass and composition as listed in table[T] The filled circles show the results at 100 Myr and 
the open triangles show the results at 10 Myr as listed in tables|2]and[3] Kepler-1 lb is extremely vulnerable to H/He mass loss, and would have been up to ~ 90% 
H/He if it formed as a water-poor super-Earth. Kepler- 11c is the least vulnerable due to its larger core mass, and d, e, and f are intermediate. All five planets are 
consistent with having initially been water-rich sub-Neptunes with comparable amounts of rock and water and ~20% of their mass in H/He. 
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Fig. 5. — Bolometric flux at the top of the atmosphere, relative to the flux incident on Earth, vs. planet density. Once again, Kepler-11 planets are shown by 
filled circles with letters indicating each planet. Open squares show the other extrasolar planets included in Figure [T] Colors indicate possible compositions. 
Planets that could be rocky are red, those that could be water-worlds are blue, and those that must have H/He are orange. For comparison, the gray crosses 
show all other transiting planets with measured masses greater than 15 and less than 100 Mq. The dashed black lines show curves of constant mass loss 
for different mass loss rates and ages, assuming our standai'd mass loss efficiency of 20 %. Finally, the shaded letters at the left indicate the densities for each 
Kepler-1 1 planet at 100 Myr predicted by our mass loss evolution models in Section [X2l 



17 



10000 F 



u 



1000 - 



100 - 



10 r 



— 1 1 — I — I — r- 



1 1 1 1 — I — I — I — r-| 

Kcplcr-u Needs H/Hr 

T f^M ' Y ''Zm Possibly Water 

+ 15 Mj, < Mp < 100 Me 

Possibly Rocky 



T 1 1 1 1 1 — I — 1— 



-m-< 

i 



— 4ffl--' 
+ + . ^ 



tioss= 12 Gyr^ 



+ 

I rn I 

J + 



+ 
+ 



_i I I I i_ 



_i I I I I I I i_ 



_i I I I I I I i_ 



1 10 100 

Mass * Density (M^ g cm"^) 



Fig. 6. — Similar to Figure |5]except liere we liave multiplied the x-axis by planet mass. Once again, Kepler-1 1 planets are shown by filled circles, while open 
squares show the other extrasolar planets included in Figure[T] Colors indicate possible compositions. Low mass planets that could be rocky are red, those that 
could be water-worlds are blue, and those that must have H/He are orange. For comparison, the gray crosses show all other transiting planets with measured 
masses greater than 15 and less than 100 Mg). There is a threshold in this diagram above which there are no observed transiting planets. Moreover, this 
threshold con'esponds to a critica l ma ss-loss timescale (see eq. |6), as shown by the dashed black line. We discuss this threshold in the context of XUV dri ven 
mass loss in section|4] In section |4~T] and figure|2]we reproduce this threshold using our coupled thermal evolution and mass loss models. Finally in section |4!2l 
and table|4]we use this threshold to provide minimum masses for the broader population of Kepler candidates without measured densities. 
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Fig. 7. — This shows the results of ^ 1000 thermal evolution and mass loss model runs which reproduce the mass loss threshold seen in Figure|6] Each panel 
plots incident bolometric flux in vs. planet density X mass for different ages and mass loss histories. Each point is sized according to its mass and colored 
according to its composition, assuming a H/He atmosphere atop an Earth-like core. The top left panel shows the initial distribution of the models before any 
mass loss has taken place. The other two top panels show results at 100 Myr and 10 Gyr for our standard mass loss efficiency (e = 0.1). Meanwhile the bottom 
panels show the results at 1 Gyr for three different mass loss efficiencies ranging from extremely inefficient (e = 0.01), to extremely efficient mass loss (e = 1). 
As planets lose mass, the points shrink, move to the right, and become bluer The dashed fine in panels 2-6, is the same as the black dashed line in Figure |6] 
corresponding to critical mass loss timesc ale. T he threshold in Figure|5]is well reproduced by models with e = 0.1, which is also the approximate value implied 
by detailed models as discussed in Section l23] 
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Fig. 8. — Separation between adjacent pairs of planets in Kepler-11, in terms of their mutual Hill spheres A. The x-axis shows the separations in terms of the 
current Hill spheres, while the y-axis shows the predicted Hill spheres when the system was 10 Myr old assuming a water-poor super-Earth scenario and that the 
planets remained stationary. The dashed lines show the approximate A > 9 stability threshold for five planet systems from Smith & Lissauer ( 2009). Likewise, 
the dotted lines show the A > 2\/3 stability threshold for two planet systems from lGladmanI fl99d/t) . Although the d-e pair does drop below the A > 9 threshold 
at 10 Myr, both the c-d and e-f pairs remain above it potentially stabilizing the system. More importantly, the b-c pair drops dangerously close to the A > 2\/3 
critical threshold for dynamical stability in two-planet systems. This is another reason we disfavor a water-poor super-Earth scenario for Kepler-1 lb. 



